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We investigate two theoretical pseudomagnon-based models for a bilayer quantum Hall system 
(BQHS) at total filling factor v t — 1. We find a unifying framework which elucidates the different 
approximations that are made. We also consider the effect of an in-plane magnetic field in BQHSs 
at vt = 1, by deriving an equation for the ground state energy from the underlying microscopic 
physics. Although this equation is derived for small in-plane fields, its predictions agree with recent 
experimental findings at stronger in-plane fields, for low electron densities. We also take into account 
finite-temperature effects by means of a renormalisation group analysis, and find that they are small 
at the temperatures that were investigated experimentally. 
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I. INTRODUCTION 

Over the last two decades, numerous experiments 
have been performed on bilayer quantum Hall systems 
(BQHSs) at total filling factor v t = 1, revealing a very 
rich physical system. In a series of groundbreaking exper- 
iments in the early years of the current decade, evid ence 
was found for the existence of an excitonic superfluidj 1 * 2 * 3 ! 
as well as of the associated Goldstone modePThis phase 
is destroyed by increasing the effective interlayer sepa- 
ration d/l, where d is the distance between the layers 
and I the magnetic length. For small effective separa- 
tions, the system is in an incompressible phase and the 
Hall conductivity exhibits a plateau at v = 1, whereas 
at large effective separations, the layers decouple and 
the Fermi-liquid behaviour of two independent layers is 
recovered; here, the system is in a compressible phase. 
The critical effective interlayer separation (d/l) c has been 
shown to be sensitive to charge imbalance, tunnelling 
amplitudes and in-plane magnetic fields. ^ In addition 
to the compressible-incompressible transition discussed 
above, a commensurate-incommensurate transition has 
been identified, driven by an in-plane magnetic field 
Recently, detailed measurements on this aspect of the 
systerrPhave become available. 

This wealth of experimental findings has naturally re- 
newed the theoretical interest in the system. At present, 
three independent models of the neutral (spin-flip) ex- 
citations of a bilayer system at Vt = 1 and the case 
of equal electron populations in each layer exist in the 
literature] 9 * 10 * 11 ! with three different derivations and two 
different predicted spectra. This state of affairs mer- 
its an investigation. That is the first aim of this pa- 
per: to compare and contrast, and where possible to 
link, the existing models. Random phase approxima- 
tion (RPA) calculations have revealed a linear Goldstone 
modep in qualitative and even rough quantitative agree- 
ment with experiment P Subsequently, an approximation 
based on pseudospin waves has been proposed for the sys- 
tem in the presence of strong tunnelling,^ which repro- 
duces the RPA result from Ref. 9 in the non-tunnelling 



limit. Lastly, in recent years, a bosonisation method has 
been proposed^ to directly study the Bose-Einstei n con- 
densate (BEC) of excitons detected experimentally. 1 234 
The same vt — 1 system with unequal layer electron pop- 
ulations has also been studied using similar techniques in 
Ref. [12l We will develop a unifying framework that links 
all three models for the balanced case and allows for fur- 
ther approximations. 

This first aim is subsidiary to a second purpose: to cre- 
ate a model for a bilayer system at v t — 1 with an in-plane 
magnetic field Bp. Apart from the above-mentioned ex- 
perimental findings, in-plane magnetic fields are often 
used to suppress tunnelling between the layers, in or- 
der to study the quantum Hall effect (QHE) in bilayer 
systems without interlayer tunnelling ]SU31 

Early theoretical work in this area has led to the 
prediction of a commensurate-incommensurate phase 
transitionPES We will revisit the problem, and derive an 
equation for the ground state energy from the underly- 
ing microscopic physics. This equation was proposed in 
Ref. |7j but is explicitly derived here. Following the work 
by Hanna et aZP^we calculate at which in-plane magnetic 
field strength the commensurate-incommensurate transi- 
tion should occur, with the aim of comparing the the- 
oretical prediction to recent experimental observations.^ 
Although this model is based on small in-plane fields, it 
turns out to agree with experiments in the regime of low 
electron densities. Including finite-temperature effects 
by means of a renormalisation group analysis (details of 
which are presented elsewhere^) yields a small change 
in critical in-plane field, which is not enough to explain 
the difference between the theoretical predictions and the 
experimental observations at larger electron densities. 

This paper is structured as follows. In section II, we 
introduce the BQHS and derive its microscopic Hamilto- 
nian. In section III, we discuss the various models that 
exist in the literature, introduce one new bosonisation 
approach, and analyse the differences between the mod- 
els. This analysis results in a unifying framework where 
all the models that are investigated here can be seen to 
be variations of each other. In section IV, we introduce 
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an in-plane field into the bilayer system, and microscop- 
ically derive an equation for the ground state energy of 
the bilayer system with in-plane field. We then present 
some results coming out of that model, including a finite- 
temperature result obtained by means of a rcnormalisa- 
tion group analysis. Section V contains conclusions and 
a discussion. 



II. THE SYSTEM 

A BQHS consists of two individual two-dimensional 
electron gas layers, parallel and at a distance d to each 
other. The QHE occurs in such systems at low temper- 
atures and under strong perpendicular magnetic fields. 
In these circumstances, the system is characterised by 
four parameters: the total filling factor i/ t , the interlayer 
separation d, the tunnelling amplitude t, and the charge 
imbalance Av. The total filling factor is given by the 
sum of the filling factors of the upper (y u ) and lower 
(vi) layers, v t = v u + Ui, whereas the charge imbalance 
Av = v u — V\ is the difference between the individual 
layer filling factors. The filling factor of an individual 
layer v — n e /n r f > counts the number of filled Landau levels 
within that layer. We will concentrate on a system with 
v t = 1 and zero charge imbalance, such that each layer 
has filling factor 1/2. This restricts the dynamics to the 
lowest Landau level (LLL), which means that we must 
project all operators into the LLL. The other two im- 
portant parameters are the effective interlayer separation 
d/l, where I — \fhfe~B is the magnetic length, and the 
ratio between the tunnelling amplitude t and the charac- 
teristic Coulomb interaction energy Eq = e 2 /el. These 
two dimensionless parameters can be tuned experimen- 
tally through the magnetic field strength B. By varying 
d/l, the relative importance of the interlayer Coulomb in- 
teraction is changed, since I is proportional to the average 
distance between two neighbouring electrons within one 
layer. Increasing t/ Eq enhances the interlayer coherence, 
since the tunnelling energy favours a coherent state. 

We will take the temperature and Zeeman splitting to 
be such that the electron spins are completely frozen. 
From a single-particle viewpoint, the layer degree of free- 
dom gives rise to a two-state system, so the relevant 
Hubert space is the same as that of a spin-1/2 system. 
Using this similarity, we map the system to a single-layer 
system where the layer degree of freedom is described by 
a pseudospin variable. In the absence of tunnelling or a 
bias voltage, there is no energy associated with the pseu- 
dospin of an electron. The Coulomb interaction, however, 
is pseudospin-dependent, since the interlayer repulsion is 
weaker than the intralayer repulsion. 

Orienting the axes of the spirP^ system such that t 
(J.) indicates an electron in a(n) (anti)symmetric super- 
position of being in both layers, the Coulomb part of the 



Hamiltonian takes the form 

H c = ^ ^2 u ^'( k )Ax,-kAx\k- (l) 

a, a' ,k 

where v a , <j'(k) is the spin-dependent Coulomb interac- 
tion. The density operator p^ ^ is given by 

PaM = e H ' k|2/2 Yl G m ,„(k)4 m c CT ,„, (2) 

where the operator c^L destroys (creates) an electron 
with spin a in the guiding centre m. The function G„ l; „ 
is defined in Appendix [X] Similarly, we can define spin 
density operators for later use: 

e -|ik| 2 /2 

S Z = 2 ^ G m,«( k ) ( C T,m C T,n ~ C l,m C l,«) > ( 3a ) 

m,n 

— |Zk| 2 /2 

Sk = — 2 £G m , n (k)(c{ im c iin + c\ m c hn ). (3b) 

The density operators defined in Eqs. ^ and (|3| obey the 
LLL projection algebra, which is discussed in Appendix 
[A"} The Coulomb Hamiltonian can be split up into total 
density and spin parts by writing 

ffc^£ «o(k)p-kPk + 2 £ w c (k)5£ k S£, (4) 

k k 

where 

«o/e(k)=^(l±e-l k l d ) (5) 

(note that we are working in a unit area system: A = 1). 
The u c -term measures the capacitive energy due to charge 
imbalance between the layers, since S x is proportional to 
the charge imbalance. Interlayer tunnelling adds a term 
Ht'' 

p e -\r\ 2 /2l 2 . 

H T = ~t d 2 r^ 2jt12 Gm,n(r)(4, m Q,n + 4,m C u,n 
m,n 

= - tS£ =0 , (6) 

where the relation between c u /i and c-fj is given by Cf / ^ — 
2~ 1 / 2 (c M ± q). Since the kinetic term is constant and 
therefore irrelevant in the v t — 1 case, we now have the 
total Hamiltonian 

H = - tS ° + «o(k)p-kPk + 2 J2 «c(k)^ k S£. (7) 

k k 

The Vo-term is invariant under rotations of the spin, 
whereas the u c -term favours S x — due to the positiv- 
ity of v c . In other words, it disfavours spin orientations 
that do not lie in the plane parallel to the layers. In the 
absence of the tunnelling term, the spin orientation is 
constrained to lie in the plane but is otherwise free. 
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Note that the orientation of the spin axes can be chosen 
freely: for example, one can choose to use spin up (down) 
to represent an electron in the upper (lower) layer. In 
that case, the roles of S x and S z in Eq. ^ are inverted 
and it is the spin z-axis that is perpendicular to the lay- 
ers. 



mllx) can be seen to be exact orthonormal eigenstates 
of the v c = Hamiltonian. Hence, our choice of spin 
orientation is appropriate to study excitations above the 
small-v c limit of the system. 



B. Models 



III. PSEUDOSPIN MODELS 



Single-mode approximation 



To accomodate the presence of a strong tunnelling 
term, we will work in the symmetric/ antisymmetric 
(S/AS) basis. Starting from a symmetric ground state, 
we study the antisymmetric excitations, which are ap- 
proximately bosonic in nature. Defining the bosonic 
vacuum to be the ferromagnet, we describe the excita- 
tions above the ferromagnet as a system of noninteract- 
ing bosons. 



A. Magnons 

In the presence of tunnelling, the symmetric state has 
the lowest energy, since the single-electron wavefunction 
for that state has no nodes in the direction perpendicu- 
lar to the layers, whereas the one for the antisymmetric 
state has a single node. The level splitting between the 
symmetric and antisymmetric states is proportional to 
the tunnelling amplitude and allows us to represent the 
ground state as a ferromagnet, in which all spins point 
in the same direction. This implies that a particular di- 
rection in spin space is selected. It is the tunnelling term 
that determines the preferred direction and thus breaks 
the in-plane symmetry. In our chosen spin orientation, 
its momentum space form is Ht = —ISq, which favours 
a uniform state in which the spin is oriented in the spin z 
direction everywhere, parallel to the plane of the layers. 

We have chosen the spin orientation such that the re- 
sulting Hamiltonian is diagonal when v c = 0, which oc- 
curs when the interlaycr distance vanishes. In that case, 
the ground state, which we denote by \x), is the state 
in which all spins are oriented along the spin z-axis: 
|x) = ||t ... t). If, on the other hand, v c ^ 0, the 
Hamiltonian is no longer diagonal, because \x) is not an 
eigenstate of S x . 

The excitations above this ferromagnet are magnons, 
created by the operator S k = S k — iS^. With the help 
of the LLL projection algebra (Appendix [A]) , it is easy 
to check that (xlKT > J S_p]|x) ^ <^p,q- This implies that 
near the ground state, modes created by S~ are approxi- 
mately bosonic. Following earlier workpnm we can define 
a magnon operator m) by normalising S~ , 



m 



D \lU\ 2 /4 



-S k , 



(8) 



where N is the total number of electrons in the system. 
By explicit calculation, the single-magnon modes |k) = 



MacDonald et al. have proposed a model based on 
the magnon-like excitations mentioned above Here, we 
briefly repeat the main results of that study with the aim 
of comparing them to the outcomes of other studies, as 
well as our own findings. The calculation of the excita- 
tion spectrum is based on a tried and tested method^ 
assuming the magnon density is low enough for the in- 
teraction to be negligible, one takes 



e k = (k|H|k) - (x\H\ X ). 



(9) 



Ignoring the magnon-magnon interaction comes down to 
computing the excitation spectrum as if there is only 
one mode; hence, it is called the single-mode approxi- 
mation. Taking into account the off-diagonal elements of 
the Hamiltonian in the magnon basis, one obtains 



H = E n 



in 



sma t m 



(10) 



\ sma 



(m k m_ k + ra k m_ k ) 



where 



= t + Ne 



-|*k| 2 /2 



w c (k) 



(11) 



-|iq| 2 /2 



A k ma =jVe-l' k l 2 /V(k) 



« c (q) + 2u (q) sin 2 (k A q/2)J , 

(12) 

w c (q)e-l Zq l 2 / 2 cos(kAq), 



£ 

q 



and we have used the shorthand notation k A q = l 2 z ■ 
(k x q) . Diagonalising Eq. ( 10 1 by means of a Bogolyubov 
transformation one obtains 



H 



(13) 



where a k are the quasiparticle excitations. The Bo- 
golyubov spectrum has the familiar form 



sma ^2 _ ^smaj2 



(14) 



0, f^ ma reduces to the spec- 



Taking the limit t 

trum found within the RPA formalism. 9 Comparing the 
above calculation of the spectrum to the one performed 
in Ref. |5J one sees that up to a minus sign, the same 
pseudospin mapping is used. The only difference is the 
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manner of calculation: the former uses Feynman dia- 
grams, whereas the latter uses the approximately bosonic 
nature of the pseudospin waves. The predicted spec- 
trum features a Goldstone mode at t — 0, as found 
experimentally^ but also a phase transition at (d/l) c rs 
1.2, whereas experiments reveal the critical separation to 
lie at (d/l) c *=s 1.8. Furthermore, the theoretically pre- 
dicted phase transition is induced by a roton minimum 
touching the axis, whereas no roton has ever been ob- 
served, in spite of an extensive search.^ 



2. S/AS bosonisation 

Recently, another approach has been presented, called 
bosonisation, which is also based on a pseudospin map- 
ping and bosonic excitations.^ However, a different spin 
orientation is used: spin up (down) refers to an elec- 
tron being in the upper (lower) layer instead of an 
(anti)symmetric state. This spin mapping is useful, since 
it allows one to describe very directly the BEC of exci- 
tons that has been observed experiment allyP^MI How- 
ever, after the bosonisation procedure, the Hamiltonian 
describes a system in which every second electron is part 
of a boson. This is problematic, because the bosonic 
operators were defined under the assumption that there 
would be few bosons: so few that they would locally only 
see the (ferromagnetic) ground state. 

Here, we will present a bosonisation scheme based on 
the S/AS splitting discussed above, but still within the 
spin mapping used in Ref . 1111 In this way, we end up with 
a system containing few bosons, maintaining the validity 
of the approximation, but the bosons being described are 
no longer the excitons observed in Refs. 1 1121.31 II 

The ground state \x) is given by J7 m c \ m \^)i an( i an_ 



tisymmetric excitations are created by the operator 



' /2 y^G m:n (k)cj m c Ti „ 



(15) 



Normalising R in order to define a proper bosonic op- 
erator, we find 



b{= e -—R-. 



N 



(16) 



Note that i? k is identical to 5 k . Now, following Ref. |2"01 
we find bosonic expressions for the operators p and S x ^ z , 
which appear in the fermionic Hamiltonian Eq. To 
obtain these expressions, we first calculate the commu- 
tators of the operators p and S z with &t. We find 



[ Pk ,61 4 ]=2z e -l (k l 2 /4 s i n ( kAq /2)6^ +k , 
[^,^]= e - |ik|2/4 cos(kAq/2)^ +k . 



(17) 



With the commutators from Eq. ( 17 1 and the actions of p 



expressions for p and S z by means of the method outlined 
in Ref. [201 

Pk = A^k,o + 2 i e -l' k l 2 /4 J2 sin(k A q/2)fe k+q 6 q , 



A^k,o 



-|ik| 2 /4 



]Tcos(kAq/2)6 k+q V (18) 



To obtain a bosonic expression for S x , we note that 
S x = (R + + R~)/2, and simply invert the definition of 
&t . Inserting the bosonic expressions for p and S x l z into 
the fermionic Hamiltonian yields a quadratic, but non- 
diagonal bosonic Hamiltonian 



H = 



\ bos 



(19) 



Again performing a Bogolyubov transformation to diago- 
nalise the Hamiltonian, we find a quasiparticle spectrum 



^bos = ^ (c bos )2 _ (A bo S)2; 



where 



A k os =Af e -l' k l 2 / 2 « c (k) and 
el° s = t + Ne-^/ 2 v c {k)+ 

2 e H ' q|2/2 Mq) sin 2 (k A q/2). 



(20) 

(21a) 
(21b) 



This result is similar to that of Sec. |IIIB~T] but not iden- 
tical, even though the same excitations were studied, and 
the same physical property (the approximately bosonic 
nature of the excitations) was the starting point for the 
approximation. In the following section, we explore the 
origin of this unexpected difference. 



C. Unification 

In the two approaches outlined above, the bosonic ex- 
citations being studied are the same (both are antisym- 
metric excitations above a symmetric ground state), and 
the method is the same (both are based on a pseudospin 
mapping), yet they yield different results. The origin of 
this apparent inconsistency can be found by expanding 
the Hamiltonian in n-magnon states: 



H =H E ^ | qi ,...,q m )x (22) 

n— qi,...,q„ m— 

(qi, • ■ ■ , q m |-ff|qm+i, ■ • ■ , qn)(qm+i, • • ■ , q n |- 



and S z on the ground state, we find the following bosonic 



Having assumed a low enough magnon density for the 
magnon-magnon interaction to be negligible, we need 
only take this series up to n = 2. From the truncated 
series, let us investigate the term J2 P q Ip) (pl^lq) (q|> or 
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more specifically, the i> c -term 



IV. IN-PLANE MAGNETIC FIELD 



2 £ « c (k)<p|SV?£|q> 



k,p,q 



= \J2 Vc(k)( X \b p (S+ k + SZ k )(S+ + s k )bl\ x ). 

k,p,q 

Expanding the brackets yields four terms, of which the 
one of interest is 



1 -|ik| 2 /2 

2 E ^^« c (k)<x|6p& k ^6 q |x). 

k,p,q 

Calculating the ground state expectation value 
(xl^p^k^k&qlx) by means of the LLL projection al- 
gebra (i.e. commuting 5*) yields 

1 -|/k| 2 /2 

2 E L ^^«c(k)(x|[[&p,[6k,bt]],& q ]|x) 

k,p,q 
1 _ e -|ik| 2 /2 

= 2 E — jv^«c(k)<x|[6p,6t]|x> 

k,p,q 



= J E e -' ik|2/2 W C (k), 



k,p,q 

which is precisely the difference between the two ap- 
proaches. In the bosonisation scheme, the commuta- 
tor [6k, &jj is replaced with its ground state expectation 
value, which is a number. Hence, the outer commutators 
vanish, and the term does not enter the spectrum. A 
similar explanation holds for the difference between Ajj os 
and A k raa . 



In short: taking the series in Eq. (22 1 up to n = 2, we 
recover the results found by MacDonald et aZ™ The dif- 
ference between the bosonisation scheme presented above 
and the one employed in Ref. EH arises from the fact 
that the former approximates the state of the system by 
the ground state already when evaluating the series in 
Eq. (22 1, whereas in the latter, this assumption is only 



made at the stage of the Bogolyubov transformation. 

It should be noted that the difference between the two 
approaches has nothing to do with the different spin rep- 
resentations used. They can be transformed into each 
other by rotating the spins, but such a rotation does not 
have any effect on the physical quantities that can be cal- 
culated from the models. The representations are only 
different in the way in which they describe the physics, 
but not in the physical approximation that is made. The 
difference in results between sections [ill B ll and lHIB 2l is 
purely a consequence of the different ways of calculating 
the bosonic Hamiltonian that were used. In the follow- 
ing, we will use the single-mode approximation described 
in Ref. EH 



Having established the method of choice, let us con- 
sider the effect of an additional in-plane magnetic field 
B\\ = _Bsin6> in the BQHS. Let us adopt a coordinate 
system where B|| = B\\ (0, — 1, 0). 24 The vector potential 
corresponding to B|| is A|| = (0, 0, cc) ; thus, a parti- 
cle tunnelling between the two layers picks up a space- 
dependent phase. The tunnelling term takes the fornfS] 



H T = - t J d 2 r h(r) • S(r), 
h(r) =(0,sin(Qa;),cos(Qx)), 



(23) 



Q = 27ri?j| d/ 0o is the characteristic momentum intro- 
duced by the in-plane field, and <p = h/e is the magnetic 
flux quantum. In momentum space, Ht is simply 



Hr = -\[S t 



Q + 



(24) 



where we have written Q = Qx. From Eq. (23 1, it is 



obvious that the tunnelling term favours spin configura- 
tions in which the spins align with h, and thus vary their 
orientation locally. Allowing for a space-dependent orien- 
tation </>(r) and then calculating the ground state energy 
yields the Pokrovsky-Talapov model, as shown below. 



A. Pokrovsky-Talapov model 

Let us construct a ground state with locally varying 
spin orientation, 



|x')=exp [i / d 2 r^(r)0(r) |x) 



exp ^E^-q)) lx>- 



(25) 



We will use the shorthand notation L = ^ q S^-q. The 
exponential rotates the spin at position r by an angle 
4>(v). It is easy to check that the rotated spin state |x') 
is properly normalised. 

We wish to calculate (x'\H\x')- Assuming that the 
spin orientation rotates slowly as a function of position, 
we can take q</>- q to be small. Expanding the exponen- 
tials, we obtain a power series in 4>(r). Taking this series 
up to O(q 2 2 ), we obtain for the Coulomb term 

U\HcW) = f / d 2 r|V0(r)| 2 , (26) 



where 



Ps 



l —r(kdkv E {k)e- {lk ^/ 2 {lk) 2 , (27) 



32tt 2 



as found in Ref. HU Here, ve is the interlayer interaction, 
given by vq — v c . For the details of the expansion, see 
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Appendix 
in Eq. pi I 



Bl Evaluating (x'\Ht\x')i with Ht as given 
we find after a similar calculation that 



(X'\H T \X') 



2irl 2 



d 2 r cos[(f>(r) - Q r] (28) 



(see Appendix B2 for the details). Combining Eqs. (28) 



and ( 26 1 , we obtain the total ground state energy in the 



presence of tunnelling and an in-plane field: 



E[4>] = J d 2 r{ 



Ps 

2 



|V0(r)| s 



t 



2ttP 



cos[0(r) — Q • r] j. 



(29) 



This is precisely the Pokrovsky-Talapov model. The 
Coulomb interaction gives rise to a spin stiffness term: in 
order to minimise the Coulomb energy, all spins should 
be aligned. The combination of a nonzero tunnelling am- 
plitude and an in-plane field results in a locally varying 
preferred spin orientation. In the case of vanishing tun- 
nelling, t = 0, all the spins will be parallel to each other 
and aligned in some arbitrary direction in the yz-planc. 
If, on the other hand, t ^ but there is no in-plane mag- 
netic field so that Q = 0, the degeneracy with respect 
to spin rotations is lifted and the spins will point in the 
z-spin direction. 



B. Commensurate-incommensurate transition 

If the spin stiffness is small, it costs little energy to 
have neighbouring spins with different orientations. In 
that case, the system will minimise the tunnelling energy 
by setting </>(r) = Q • r, i.e. the spins follow the rotation 
imposed by the tunnelling term. This is called the com- 
mensurate phase. If, on the other hand, the spin stiffness 
is large, the gradient term represents a high energy cost 
associated with a non-uniform spin orientation and the 
system will give up the tunnelling energy in favour of a 
better Coulomb correlation. In that case, the spins do not 
follow the tunnelling term; this is called the incommen- 
surate phase. The onset of this phase is characterised by 
the appearance of solitons: sudden rotations of the spin 
orientation by 2ir. Fig. 1 shows a sketch of the behaviour 
of (f>(r) in the incommensurate phase, with two solitons 
visible. 



cf>(x) 



FIG. 1: (Colour online) Sketch of (f>(x) in the incommensurate 
phase. The gray (dashed) trace indicates Qx while the blue 
(solid) indicates 4>(x). Two solitons (sudden rotations of the 
phase by 2n) are visible in the figure. In the commensurate 
phase, <f) would exactly follow Qx. 



rotation by 2ir. The phase transition occurs when a fi- 
nite soliton density is energetically favourable compared 
to the commensurate (zero- solito n) state. The critical 
value of \IQ\ is then given by! 15 * 16 1 



|JQc| = 



3/2 



(30) 



Recalling that -By = B±t&n9 and that tan# = \lQ\l/d, 
we find 



B«=B\ 



I 

d \tt 



3/2 



t 



(31) 



Since p s depends on d/l and Eq, which in turn depend 
only on B± and constants, this equation predicts the crit- 
ical in-plane field for a given perpendicular field. In or- 
der to compare the theory to the experimental findings 
presented in Ref. [HI we need to express the total electron 
density ut as a function of the critical in-plane field. The 
total electron density is related to B± by 



1 



2ttZ 2 2ttH ! 



(32) 



so we can simply invert Eq. (31 1 to obtain nx as a func 
tion of Bf, for the comparison. 



1. Zero temperature 



2. Finite temperature 



The commensurate-incommensurate transition is gov- 
erned by two quantities: the ratio t/p s and the modulus 
of the characteristic momentum |ZQ| = (d/l) tan 6. In 
order to compare the theory to experiments, we need an 
equation for the critical in-plane field strength. Neglect- 
ing finite-temperature effects, we obtain such an equation 
from the energy functional given in Eq. ( 29 1 . We consider 



To include the effect of a finite temperature, one has 
to consider the partition function^ 



V<t>e 



-E[<t>]/k B T 



(33) 



where fc^ is Boltzmann's constant and T the tempera- 
ture. The functional integral in Eq. (33 1 can be per- 



the energy of a single soliton in the system, i.e. a single formed step by step, using the renormalisation group 



(RG) technique. In every step, an infinitesimal part of 
the integration is carried out, yielding an effective energy 
functional at an infinitcsimally lower momentum scale, or 
cut-off. By computing the change in the parameters in 
the energy functional under an infinitesimal change in 
the cut-off, one obtains the flow equations for p s and tP^ 
Integrating these equations over the cut-off running from 
its initial value to zero, one obtains an effective energy 
functional E'[4>] for which we have 



-E'[4,]/k B T 



(34) 



E' contains the effective values of p s and t at temperature 
T, allowing us to compute the renormaliscd value of the 
critical |ZQ| 



\lQc 



2\ 3/2 

Ps(0) ' 



(35) 



where e = ln(Ao/A), with Ao being the initial value of the 
cut-off A. The RG procedure for the Pokrovsky-Talapov 
model is discussed in detail in Ref. [TBI 
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FIG. 2: (Colour online) In-plane field vs. total electron density 
phase diagram indicating the commensurate-incommensurate 
phase transition. The experimental data are given by the 
red dots (the trace is a guide for the eye). The theoretical 
predictions are given by the blue, dashed line (mean field) and 
the green, solid line (finite temperature). The experimental 
data are taken from Ref. [8] 



C. Comparison with experiment 



In recent experimental work, a commensurate- 
incommensurate phase transition has been accurately 
measuredP Working at a temperature T = 130 mK, mea- 
suring on a sample with d = 23 nm and t = Asas/2 = 5.5 
K, the red dots in Fig. 2 were obtained. The parameter 
values reported in Ref. .3 give a |ZQ C | of about 2.1-2.2. 
This value invalidates the assumption that |ZQ| -C 1, 
which was made to derive the Pokrovsky-Talapov model. 
Nonetheless, it turns out that Eq. (31 1 predicts the exper- 



imentally observed critical in-plane fields in the regime 
of low electron densities, as can be seen in Fig. 2. As 
the electron density increases, the accuracy of the pre- 
diction decreases. This may be related to the fact that 
increasing the electron density is equivalent to decreasing 
the magnetic length I, and hence, to increasing the effec- 
tive interlayer separation d/l, whereas the dependence of 
p s on d/l was derived under the assumption that d/l is 
small. 

Including finite-temperature effects in the manner de- 
scribed in Sec. |IV B 2| reduces the predicted critical in- 
plane field (see Ref. |T5J) . At the temperature reported 
in Ref. 2 we find a very small change in the critical in- 
plane field strength (see Fig. 2, green trace) . On the other 
hand, the smallness of the shift validates the assumptions 
made during the experimental data analysis.^ 

Testing the theory on the experimental data provided 
by Murphy et alW, who measured on a sample with ny = 
1.26 x 10 15 m~ 2 , t = 0.4 K and d = 21 nm, we find a 
critical angle of inclination of 18°, about a factor of 2 
off from the reported value of 8 ± 2°, as already found 
by Kun Yang et aZ.P^ Taking finite-temperature (T = 0.4 
K) effects into account, this prediction improves to 15.7°: 
again a small change in the right direction. 



DISCUSSION & CONCLUSIONS 



In this paper, our aim was twofold: firstly, to clarify 
the situation of the many magnon models for a BQHS at 
i>t = 1; and secondly, to study the effects of an in- plane 
field in such a system. We have considered two models for 
a BQHS at zero temperature and vt = 1, both of which 
turn out to be variations of the magnon model obtained 
by expanding the Hamiltonian in n-magnon states. We 
have also seen what the origin of the difference between 
the two models is: it is the moment in the calculation at 
which the magnons are approximated to be bosonic. 

We have derived the Pokrovsky-Talapov model for the 
ground state energy of a bilayer QHS at Vt = 1 with an 
in-plane field from the underlying microscopic physics, 
which had been suggested earlier by Kun Yang et alV^ We 
have seen the phases predicted by this model in the weak 
in-plane field limit: a commensurate phase, in which the 
pseudospin follows the underlying structure provided by 
the combination of the tunnelling term and in-plane field, 
and an incommensurate phase, where the pseudospin ro- 
tates incommcnsuratcly with the underlying structure. 
We observed that this model, although derived for weak 
in-plane fields, predicts experimentally measured values 
in the regime of low electron densities, and we have of- 
fered an explanation for the loss of agreement at higher 
electron densities. We have showed that the inclusion of 
finite-temperature effects produces only small effects at 
the temperatures reported in Refs. 16181 Even though the 
effects considered here are small, a finitc-tcmperature- 
adjusted theory is a potentially valuable asset in future 
research. 
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APPENDIX A: THE LLL PROJECTION 
ALGEBRA 

The function G m „(k) is given by^ 



G m ,„(k) = 9(m - n)<\ 



-ilk* 



T ' 



1*1 



+ 9(n — m)\ 



,1 



—ilk 



T '' 



|/k| 



(Al) 



where k — k x + ik y . Making use of the Landau level 
basis, we can calculate the sum of the product of two 
G m ,„-functions (i.e., the matrix product G(k)G(q)). We 
obtain 

G„u(k)G,„(q) = e- i2k ^ 2 e- lkA ^ 2 G„ v(1 (k + q) 

(A2) 



where k A q = I 2 z ■ (k x q) . Eq. ( A2 1 gives rise to the LLL 
projection algebra, which can be summarised as follows: 

[Pk, P«J = e z2k ' q/2 2i sin(k A q/2)p k+q 
[^^q]=e' 2k - q/2 2 i sin(kAq/2)^ +q 

[S£, S$ = e^ 2 sin(k A q/2)p k+q (A3) 

+ ie'" /CT cos(kAq/2)S k r +q 

We can also derive the following commutators, for later 
use: 

K ,S- } = e fu ^ 2 (i sin(k A q/2)p k+q 

+ 2cos(kAq/2)S k : +q 



[Sk.Ai] -e' 2k 'i/ 2 2isin(kAq/2) > S i 
[S£,S'\ = ± e 1 ^ 2 cos(k A q/2)S ; 



(A4) 



k+q 
k+q- 



APPENDIX B: EXPECTATION VALUES IN THE 
COMMENSURATE (ROTATING-SPIN) GROUND 
STATE 

We are interested in expectation values of the form 
(x'|Ak|x')j where \x') is the rotating-spin ground state 



defined in Sec. |IV A| Hence, we need to evaluate terms 
of the form e~ lL A-^e lT . To this end, we use a simple 
generalization of the Baker-Hausdorff lemma, and find 



-iv 



A k e ir = ^G n (^ k ) 



(Bl) 



where 



■ qi,...,q» (B2) 



Coulomb term 



Using Eq. (A3 1, it is easy to verify that for odd n, 
C n (pk) oc 5g +qi+ ... +qn , whereas for even n, G„(p k ) cx 

p k+qi ^ |_q n . For C n (S^.) 1 we obtain the same solutions, 

but even and odd n are inverted with respect to G„(p k ). 
Hence, projecting the Coulomb term He into the spin- 
rotating ground state \x'), we obtain 

+ 2 J2 v c (k)(x'\C n (S%)C nl (SZ)\x'). 



In Ref. EU this series is taken up to n + n' = 2. In 
that case, with the assumption of a slowly varying spin 
orientation, one obtains the result presented in Eq. (26 1. 



9 



2. Tunnelling term 

The tunnelling term Ht can be analysed in the same way. We need 

(x'\h t \x') = -\ £<xo| (c^) + c n (s*_ Q ) + ic n (sv Q ) - iC n (sy_ Q )) \ X o). 



(B3) 



From E q. flA3| , we can deduce that the surviving terms are of the form C2n(S^) and i C2n+i(S^.) ■ Repeated application 
of Eq. QA3P yields 



(-iy 

(2n)\ 



E 



J 2 q 2n -Q/2 J 2 qr(q2 + --+Q)/2 



(B4) 



x cos(q 2 „ A Q/2) . . . cos(qi A (q 2 + ■ • • + Q)/2)S , q +qi _ 



and similarly for C^m+iXSq). Since we are assuming a slow modulation, we are working at small ZQ, and we may take 

Z 2 |Q||q| w 0. Evaluating Sq +qi _| ^_ t . n in the uniform ground state gives a Kronecker delta; inserting this into the 

sum from Eq. (B4| and Fourier transforming back to real space, we find 

J n [ >■ 

The sum is an expansion in powers of <f> of the function cos(</>(r)). Similar calculations for the other three terms from 
Eq. (pi yield 



(M^M) =- N 2j ^{e lQ ' r [cos(0(r)) - isirn>(r)) 
= - J d 2 r cos(0(r) - Q r). 



,-*Q-r 



cos(0(r)) + isin(0(r)) | 



(B5) 



This equation appeared first in Ref. 13 but since it was not derived there, we include it here. 
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Since the real electron spins do not play a role in the sub- 
sequent discussion, we will drop the prefix 'pseudo' and 
refer to pseudospin as 'spin'. 

This somewhat strange choice of coordinate system is made 



to facilitate the comparison with Ref. [14] Note that from 
this point on, we are again using the spin orientation cho- 
sen in Sec. [II] 



